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ABSTRACT 

Passive early-type galaxies (ETGs) provide an ideal laboratory for studying the interplay between dust formation around evolved stars and its 
subsequent destruction in a hot gas. Using Spitzer-IRS and Herschel data we compare the dust production rate in the envelopes of evolved AGB 
stars with a constraint on the total dust mass. Early-type galaxies which appear to be truly passively evolving are not detected by Herschel. We 
thus derive a distance independent upper limit to the dust grain survival time in the hostile environment of ETGs of < 46 ± 25 Myr for amorphous 
silicate grains. This implies that ETGs which are detected at far-infrared wavelengths have acquired a cool dusty medium via interaction. Given 
likely time-scales for ram-pressure stripping, this also implies that only galaxies with dust in a cool (atomic) medium can release dust into the 
intra-cluster medium. 

Key words. Galaxies: elliptical and lenticular, cD - Galaxies: ISM - Infrared: galaxies 
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1. Introduction 

In this paper we compare the dust outflow rate from the evolved 
AGB stars in passive ETGs, with the total dust mass, to esti- 
mate dust grain lifetimes. Both quantities are derived from ob- 
servations; the outflow rate from Spitzer-IRS spectra, and the to- 
tal dust mass from Herschel (Pilbratt et al. 2010) maps of the 
Herschel Virgo Cluster Survey, HeViCS (Davies et al. 2010, 
www.hevics.org). We assume that evolved stars are the only 
source of dust in these passive systems as there is very little evi- 
dence that type la supernovae produce dust. 

IRAS data have been used in the past in a similar way. Soifer 
et al. (1986) estimated the mass-loss from evolved stars in the 
bulge of M 31 from 12 and 25 pm fluxes and the total dust mass 
from the 60 and 100 pm fluxes. They concluded that given that 
the observed dust mass could be produced by stellar outflows in 
only 10 7 yr, there must be some mechanism that depletes the 
bulge of inter-stellar matter. Jura et al. (1987) made similar cal- 
culations for a sample of elliptical galaxies. They estimated dust 
outflow rates by assuming that half of the 12 pm flux came from 
dusty stellar envelopes. They also found that the stellar mass- 
loss rate would be sufficient to produce the observed cool ISM 
in much less than a Hubble time. 

Dust mass production rate and total dust mass can only rea- 
sonably be compared in objects that are truly passively evolv- 
ing. Temi et al. (2007) find, that for a given blue luminosity, 
the 70 and 160 pm luminosities vary by 2 orders of magnitude, 
indicating that many ETGs have a significant dust component 
not directly attributable to mass-loss from evolved stars. For our 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



purposes, these ETGs are not 'passive' because other processes 
(mergers?) have probably contributed to the dust mass. 

Recent studies of ETGs with Spitzer (Bressan et al. 2006; 
Clemens et al. 2009) have shown that even in samples of ETGs 
selected to be the most passive objects (ie. lying on the colour- 
magnitude relation) a significant fraction show evidence of either 
on-going or recent past star formation. These results indicate that 
ETGs in which dust features are seen in the optical (e.g. Sadler 
& Gerhard, 1985) and many that have been detected by IRAS at 
60 and 100 pm (Knapp et al. 1989) and Spitzer at 70 and 160 pm 
(Kaneda et al. 2007) may also host low levels of star formation. 
Intriguigly, of the 7 ellipticals detected at 70 and 160 pm by 
Kaneda et al. (2007) the only object not showing PAH emission 
features is actually a radio source, so that the far-infrared emis- 
sion is probably synchrotron! 

Tsai and Mathews (1995) studied dust destruction in ETGs 
via thermal sputtering in the hot (10 6 - 10 7 K), low density 
(n{/ -0.1 cm 4 ) gas, and found that the destruction time-scale is 
short compared to any cooling flow or dust transport time-scale. 
They therefore concluded that dust is destroyed 'on the spot' be- 
fore it has time to migrate within a galaxy. In a following work, 
Tsai and Mathews (1996) found dust-to-gas mass ratios that are 
orders of magnitude less than that in typical stellar ejecta or in 
the ISM of the Milky Way. Their interpretation of the IRAS data 
was, at that time, hampered by the resolution of the existing ob- 
servations; they were therefore not able to tie down the origin 
and distribution of the observed dust emission. 

As part of the Herschel Science Demonstration Phase (SDP), 
a 4 x 4 deg 2 field, centred approximately on M 87, has been 
observed as part of the HeViCS with both the PACS (Poglitsch et 
al. 2010) and SPIRE (Griffin et al. 2010) instruments. See Davies 
etal. (2010) for details. 
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2. Results 

Bressan et al. (2006) observed a sample of ETGs in the Virgo 
cluster with Spitzer-IRS. The objects were selected to lie on the 
optical colour-magnitude relation and were expected to be pas- 
sively evolving objects. The mid-infrared spectra for this sample 
showed that only 76% were actually passive objects, with the 
remainder showing signs of either on-going or recent past star 
formation (PAH features) or AGN activity (atomic lines). Of this 
sample of 17 galaxies, 9 lie in the HeViCS SDP field. 6 of these 
show IRS spectra consistent with passively evolving stellar pop- 
ulations (with the possible exception of NGC 4371), and it is 
these 6 objects, listed in Table [1] on which we base our study. 

2.1. Dust mass detection limit 

Of the 6 passive ETGs in the HeViCS SDP field, none are de- 
tected at any wavelength (100, 160, 250, 350 and 500 //mfl 
Because of the sensitivity of the emission at the shorter PACS 
wavelengths to small masses of relatively warm dust, we con- 
sider only the SPIRE images at 250, 350 and 500 pm to deter- 
mine the dust mass detection limit. 

Our mass detection limit depends on whether we consider 
the detection of an extended or point source. For passive galax- 
ies, the dust should follow the stellar distribution, and thus be 
strongly centrally peaked. As the core radius of an elliptical 
galaxy at the distance of the Virgo cluster is typically less than 
the resolution of the SPIRE images (respectively 18", 25" and 
37") the most likely detection would take the form of a point 
source. 

Unless the dust is very cold (< 17 K) the tightest limit on 
the dust emission comes from the 250 pm map. This has a point 
source detection limit of S lms =9.6 mJy. We calculate the dust 
mass detection limit as (Hildebrand, 1983): 



M d < 



s rm s(y)D 2 

K d (v)B(v, 7d) 



(1) 



where D is the distance to the Virgo cluster (16.5 Mpc, Mei 
et al. 2007) and ka(v) is the dust mass opacity coefficient. 
The value of Kd(v) depends on the grain composition but the 
mid-infrared spectra (§ 12. 2t are consistent with pure silicate 
grains. For dust composed only of amorphous silicate grains 
* d (250 pm) = 0.6517 m 2 kg- J (Draine & Lee, 1984). Adopting 
this value, and Tj = 30 K, we arrive at a 2cr mass detection limit 
ofM d < 8.7x1 3 MeB 

2.2. Dust production rate from AGB stars 

Assuming a spherically symmetric, stationary stellar wind of ve- 
locity Voo, and dust mass-loss rate M CE the circumstellar dust 
density is, 



p7 



M CE 
Anr 2 v a 



(2) 



The optical depth of the circumstellar envelope, r, is ob- 
tained by multiplying Eq. (0 by the dust opacity coefficient and 



1 We note that several ETGs are detected in the SDP field, including 
M 87, M 84, M 86, NGC 4459, NGC 4435 and NGC 4550, but none of 
these are truly passively evolving. 

2 Were the dust composed of amorphous carbon (Serra Diaz-Cafio & 
Jones, 2008) then * d (250 pm) = 1.932 (Rouleau & Martin, 1991) and 
the mass detection limit M d < 2.9 x 10 3 M Q . 



integrating from r,„ to r out , where the former is the dust sublima- 
tion radius and the latter is a convenient outer radius >> r,„. We 
thus have, 



T\ 



M CE k x 
Anv^r;,, 



(3) 



where k\ is the dust opacity at Ipm, following the notation of 
Bressan et al. (1998), and we have neglected the very small term 
\jr out . The dust sublimation radius is r,„ oc L 1 ^ 2 (Granato & 
Danese 1994) and the proportionality constant depends mainly 
on the dust composition (through its sublimation temperature), 
the size distribution and only weakly on the shape of the stellar 
spectrum. Inserting the value of r,„into Eq. (O we obtain, 
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where L\ is the bolometric luminosity in units of 10 4 L Q . The 
factor A is mainly determined by the dust mixture. For the mix- 
ture of silicate grains of M-type AGB stars used in Bressan et 
al. (1998), suitable for the spectra of passive ETGs (Fig. [TJ, 
A^ 6.5 x lO- IO M yr I /(kms- 1 ). 

To determine the dust mass-loss rate from the MIR spec- 
trum we fit a combination of a pure photospheric atmosphere, 
a MARCS model (Gustafsson et al. 2008) of 4000 K, plus the 
emission from a dusty envelope selected from the library of 
Bressan et al. (1998), as described in Bressan et al. (2007). We 
thus obtain the fraction of flux, say at 10;um, /^(lOpm), due 
to dusty circumstellar envelopes, and, from the distance of the 
galaxy, the corresponding "dust" luminosity sampled by the slit, 
AnD 2 ff\\0pm). 

Surprisingly, though ETGs contain a mixture of evolved stars 
of different metallicity and age, a single dusty envelope is typi- 
cally enough for a good fit to the MIR spectrum. 

The shape of the broad feature near 10/im in general requires 
T\ ~ 2-A, and so the envelopes are optically thin at 10/im. From 
the shape of the 10pm feature we thus derive the optical depth 
of the "average" AGB star. To compute its dust mass-loss rate 
through Eq. (0|i we need to assume a bolometric luminosity and 
wind velocity. Typical values for evolved, old stars are La, ~0.2- 
0.4and Voo -5-15 kms -1 . 

The total dust mass-loss sampled by the slit is finally ob- 
tained by multiplying the value for the "average" AGB star by 
the number of dusty AGB stars sampled by the slit. This is given 
by the ratio of the dust luminosity of the galaxy to that of the 
"average" AGB star, L d E (10pm). Thus, the total dust mass-loss 
sampled by the slit is, 



M 



slit 



At\ Vca VL4 



AnD 2 ff\\0pm) 
LfilOpm) 



(5) 



Summarizing, the spectral fit provides t\ and f d ll (l0pm). 
From n, which characterizes the dusty envelope, and the as- 
sumed bolometric stellar luminosity, we obtain L d E (10pm) and 
finally M d lt , for a specified value of v^. We note that, since 
L d E (lQpm) scales approximately with the bolometric luminos- 
ity of the star, the total mass-loss has only a weak dependence 
on the assumed luminosity of the typical AGB star. 

We finally note that recent work (Boyer et al. 2010) has 
shown dust to be released into the ISM only by AGB stars, con- 
trary to some claims that RGB stars also contribute. 
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Fig.l. Spitzer-IRS spectrum of NGC 4473. The fitted model 
components are the stellar photosphere, and dusty circumstel- 
lar envelope (CE). The spectrum of NGC 4371 is shifted by 0.2 
dex. 

Table 1. Dust mass loss rates and grain lifetimes, f g . 



NGC 


MJL V (KY 


M i /M t 


M d (beam) 


1 




Moyr -1 
10 26 WHz-' 


Msyr -1 

10 12 Afo 


10 


4 MqVt- 1 


10 6 yr 


4371 


3.80 


0.00855 




1.0 


86 


4442 


3.48 


0.00712 




1.7 


51 


4473 


3.57 


0.00742 




1.9 


46 


4474 


3.71 


0.00806 




0.59 


150 


4551 


2.61 


0.00468 




0.36 


240 


4564 


4.28 


0.00959 




1.3 


65 



Notes. * For an assumed dust temperature of 30 K. f The specific dust 
mass-loss rate per unit K-band luminosity. * The mass loss per 10 12 M Q 
of galaxy mass. A dust mass opacity coefficient suitable for amorphous 
silicate grains has been used. 



2.3. Dust grain lifetimes 

Dust mass-loss rates were derived from Spitzer-IRS spectra that 
were taken through a slit of dimensions 18"x 3'.'6 centred on the 
galaxy. In order to compare these values with the dust masses 
within one 250 pm beam we assume that the light distribution 
in the mid and far-infrared are similar, and that both are similar 
to that at K-band. For each galaxy, we use the 2MASS K-band 
image, smoothed to the Spitzer-IRS spatial resolution, to deter- 
mine the flux within the IRS slit. We then smooth the K-band 
image to the Herschel 250 pm resolution and calculate the flux 
within one beam. We use the ratio of these two fluxes to derive 
the mass-loss rate from the area of one Herschel beam. 

We derive dust grain lifetimes simply as the ratio of dust 
mass detection limit to dust mass-loss rate, Md/Mj. Our dust 
mass detection limit of M& < 8.7 x 10 3 M Q for silicate grains 
results in grain lifetimes as given in Table Q] The tightest con- 
straint, of < 46 + 25 Myr, is provided by NGC 4473. The error 
reflects the uncertainties in v ro and L4. In Table|2]we provide the 
values used to arrive at this grain lifetime and their dependen- 
cies. For dust temperatures in the range 20 - 40 K, for example, 
the grain lifetime lies in the range 69 - 3 1 Myr. 

If dust grains are destroyed via sputtering in the 10 6 - 10 7 K 
gas in ETGs then the theoretically, expected dust grain lifetime 
in years is t — kj r g /rin (e.g., Tielens et al. 1994), where r g is the 
grain radius in nm, «h is the gas density in cirT 3 and k[ = 310 
for silicate dust (210 for amorphous carbon, 1500 for graphite) 



Table 2. Parameter values for calculation of grain lifetime and 
the power of their dependence on the derived lifetime. 



»J mis 

mjy 



*d(250) 
m 2 kg -1 



K 



A* 

10- 10 



T\ 



kms 



u 

10 4 L R 



9.6 
1 



0.6517 
-1 



30 
-1+ 



6.5 

1 



fit 
1 



10 
1 



0.3 
0.5 



Notes. E.g., t g oc l/r d . * The units of A are M Q yr '/(kms '). f In the 
Rayleigh- Jeans limit; which is a bad approximation for Ta < 20 K. 



so the grain lifetime should be longer in galaxies with a lower 
density hot ISM. 

Of our sample, NGC 4371 and NGC 4474 both lie in the 
least X-ray luminous category as defined by Irwin & Sarazin 
(1998). NGC 4564 has been observed by Chandra, and the cen- 
tral hot gas density estimated to be 0.011 + 0.006 crrT 3 by 
Soria et al. (2006). These authors analyzed a sample of 'quies- 
cent' galaxies and found central gas densities in the range 0.01 1- 
0.03 crrT 3 . These values overlap with the low end of estimates 
from Chandra data of more general ETG samples (Pellegrini, 
2005). We thus assume a hot gas density of 0.02 crrT 3 for all our 
ETGs. The expected dust grain lifetime, for "typical" 100 nm 
silicate grains is 1.55 Myr (1.05 for amorphous hydrocarbon, 
7.5 for graphite) and thus consistent with our non-detections. 



3. Discussion and conclusions 

We interpret the ratio of dust mass to dust production rate as a 
measure of grain lifetime. In principle, however, dust could be 
removed from the galaxy rather than destroyed. 

A wind of hot gas produced by type la supernovae probably 
cannot drive dust (or gas) out of the potential well of massive 
ETGs (Mathews & Loewenstein 1986). However, good evidence 
of ram-pressure stripping of the hot ISM in elliptical galaxies 
has been found in the form of X-ray tails (Randall et al. 2008; 
Machacek et al. 2006). If dust is found in this hot gas then it 
too may be removed. However, the timescale for such stripping 
is estimated to be of the order a few 10 8 yr (Takeda et al. 1984; 
Murakami & Babel 1999, eq. 12) which is an order of magnitude 
longer than our upper limit to the grain lifetime. The observa- 
tion that the X-ray emission of most ETGs is not displaced from 
the stellar distribution, in fact, suggestes that either the strip- 
ping time-scale is longer than the hot gas production time-scale 
or that the ram-pressure is typically insufficient to remove the 
hot medium. It is now clear that if dust is isolated from the hot 
medium in denser, cooler clouds, than ram-pressure does have 
time to strip dust (Cortese et al. 2010) but the passive ETGs in 
our study do not show any evidence of accumulating dust in such 
clouds. Passive ETGs do not directly pollute the intra-cluster 
medium with dust. 

Grain lifetimes of < 46 Myr have immediate implications 
for ETGs that are detected in the FIR. Either they produce dust 
at rates much higher than those given in Table Q] or their dust is 
shielded from the hot gas in cool clouds. For example, the dust 
mass of NGC 4435 of 1.2 x 10 6 M Q would require a total dust 
production rate of 0.03 M yr _1 , more than an order of magni- 
tude greater than any value in Table Q] This object is known to 
host star formation (Panuzzo et al., 2007) and so dust is likely to 
be in cool clouds more typical of late-type galaxies; grain life- 
times in this case may be much longer. Dust may have been ac- 
quired from an interaction, but in any case it probably arrived as 
part of a cooler medium. 
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In a passive ETG, the ISM, and the dust within it, come 
mainly from mass-loss from evolved stars. Unless the dust-to- 
gas ratio in this released material is very variable, one would 
expect dust and gas production to follow one another. If the 
X-ray luminosity, L x oc n^ one would naively expect a rela- 
tion L x oc M?. However, if the grain lifetime, t s oc l/n H , then 
we would expect only a linear correlation between L x and M&. 
Kaneda et al. (2007) find evidence of an anti-correlation for a 
small sample of elliptical galaxies detected by Spitzer-MlPS, but 
this may be expected if the dust were of external origin. 

Only NGC 4371 is detected by MIPS at 70 pm (marginally 
at 160 yum). The detection is consistent with our Herschel de- 
tection limits. This object actually shows extremely weak PAH 
features in its Spitzer-IRS spectrum, and is therefore probably 
not totally passive. The object is included for comparison, and 
its IRS spectrum included in Fig.[T] 

We have assumed that AGB stars are the only source of dust. 
Were type la supernovae able to produce dust from the metals 
they produce, we estimate that the contribution to the dust pro- 
duction rate would be, at most, similar to that of AGB stars. Our 
limit on grain lifetime would be a factor of ~2 lower in this case. 

Although dust grains are rather short-lived in ETGs, the 
prospects for survival in the intra-cluster medium are rather 
better because the gas densities are typically at least 2 orders 
of magnitude lower than in ETGs (Bohringer et al., 1994). 
Therefore, if dust can be removed from cluster galaxies, or form 
in the intra-cluster medium, there is a good chance that such dust 
will be detected by the HeViCS, especially at larger cluster radii 
where gas densities are lower. 

Acknowledgements. MC and AB acknowledge support from contract ASI/INAF 
1/016/07/0. 
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